Dispersal Patterns and Influence of Air Travel During the Global Expansion of SARS-CoV-2 Variants of Concern

The Alpha, Beta and Gamma SARS-CoV-2 Variants of Concern (VOCs) co-circulated globally during 2020-21, fueling waves of infections. They were displaced by Delta during a third wave worldwide in 2021, in turn displaced by Omicron in late 2021. In this study, we use phylogenetic and phylogeographic methods to reconstruct the dispersal patterns of VOCs worldwide. We find that source-sink dynamics varied substantially by VOC, and identify countries that acted as global and regional hubs of dissemination. We demonstrate a declining role of presumed origin countries of VOCs to their global dispersal, estimating that India contributed <15% of Delta exports and South Africa <1-2% of Omicron dispersal. We estimate that >80 countries had received introductions of Omicron within 100 days of emergence, associated with accelerating passenger air travel and higher transmissibility. Our study highlights the rapid dispersal of highly transmissible variants with implications for genomic surveillance along the hierarchical airline network. Graphical Data analysis clarifies that dispersal of SARS-CoV-2 variants from their sites of initial detection was related to the amount of global air travel at the time of the variant's emergence, and that travel volume through "hub” sites distinct from the site of emergence was a key driver of variant spread.

Forecasting hospital-level COVID-19 admissions using real-time mobility data.

BACKGROUND: For each of the COVID-19 pandemic waves, hospitals have had to plan for deploying surge capacity and resources to manage large but transient increases in COVID-19 admissions. While a lot of effort has gone into predicting regional trends in COVID-19 cases and hospitalizations, there are far fewer successful tools for creating accurate hospital-level forecasts. METHODS: Large-scale, anonymized mobile phone data has been shown to correlate with regional case counts during the first two waves of the pandemic (spring 2020, and fall/winter 2021). Building off this success, we developed a multi-step, recursive forecasting model to predict individual hospital admissions; this model incorporates the following data: (i) hospital-level COVID-19 admissions, (ii) statewide test positivity data, and (iii) aggregate measures of large-scale human mobility, contact patterns, and commuting volume. RESULTS: Incorporating large-scale, aggregate mobility data as exogenous variables in prediction models allows us to make hospital-specific COVID-19 admission forecasts 21 days ahead. We show this through highly accurate predictions of hospital admissions for five hospitals in Massachusetts during the first year of the COVID-19 pandemic. CONCLUSIONS: The high predictive capability of the model was achieved by combining anonymized, aggregated mobile device data about users' contact patterns, commuting volume, and mobility range with COVID hospitalizations and test-positivity data. Mobility-informed forecasting models can increase the lead-time of accurate predictions for individual hospitals, giving managers valuable time to strategize how best to allocate resources to manage forthcoming surges.

During the COVID-19 pandemic, hospitals have needed to make challenging decisions around staffing and preparedness based on estimates of the number of admissions multiple weeks ahead. Forecasting techniques using methods from machine learning have been successfully applied to predict hospital admissions statewide, but the ability to accurately predict individual hospital admissions has proved elusive. Here, we incorporate details of the movement of people obtained from mobile phone data into a model that makes accurate predictions of the number of people who will be hospitalized 21 days ahead. This model will be useful for administrators and healthcare workers to plan staffing and discharge of patients to ensure adequate capacity to deal with forthcoming hospital admissions.

Global disparities in SARS-CoV-2 genomic surveillance.

Genomic sequencing is essential to track the evolution and spread of SARS-CoV-2, optimize molecular tests, treatments, vaccines, and guide public health responses. To investigate the global SARS-CoV-2 genomic surveillance, we used sequences shared via GISAID to estimate the impact of sequencing intensity and turnaround times on variant detection in 189 countries. In the first two years of the pandemic, 78% of high-income countries sequenced >0.5% of their COVID-19 cases, while 42% of low- and middle-income countries reached that mark. Around 25% of the genomes from high income countries were submitted within 21 days, a pattern observed in 5% of the genomes from low- and middle-income countries. We found that sequencing around 0.5% of the cases, with a turnaround time <21 days, could provide a benchmark for SARS-CoV-2 genomic surveillance. Socioeconomic inequalities undermine the global pandemic preparedness, and efforts must be made to support low- and middle-income countries improve their local sequencing capacity.

Shifts in global mobility dictate the synchrony of SARS-CoV-2 epidemic waves.

BACKGROUND: Human mobility changed in unprecedented ways during the SARS-CoV-2 pandemic. In March and April 2020, when lockdowns and large travel restrictions began in most countries, global air-travel almost entirely halted (92% decrease in commercial global air travel in the months between February and April 2020). Initial recovery in global air travel started around July 2020 and subsequently nearly tripled between May and July 2021. Here, we aim to establish a preliminary link between global mobility patterns and the synchrony of SARS-CoV-2 epidemic waves across the world. METHODS: We compare epidemic peaks and human global mobility in two time periods: November 2020 to February 2021 (when just over 70 million passengers travelled), and November 2021 to February 2022 (when more than 200 million passengers travelled). We calculate the time interval during which continental epidemic peaks occurred for both of these time periods, and we calculate the pairwise correlations of epidemic waves between all pairs of countries for the same time periods. RESULTS: We find that as air travel increases at the end of 2021, epidemic peaks around the world are more synchronous with one another, both globally and regionally. Continental epidemic peaks occur globally within a 20 day interval at the end of 2021 compared to 73 days at the end of 2020, and epidemic waves globally are more correlated with one another at the end of 2021. CONCLUSIONS: This suggests that the rebound in human mobility dictates the synchrony of global and regional epidemic waves. In line with theoretical work, we show that in a more connected world, epidemic dynamics are more synchronized.

Preparedness and vulnerability of African countries against importations of COVID-19: a modelling study.

BACKGROUND: The novel coronavirus disease 2019 (COVID-19) epidemic has spread from China to 25 countries. Local cycles of transmission have already occurred in 12 countries after case importation. In Africa, Egypt has so far confirmed one case. The management and control of COVID-19 importations heavily rely on a country's health capacity. Here we evaluate the preparedness and vulnerability of African countries against their risk of importation of COVID-19. METHODS: We used data on the volume of air travel departing from airports in the infected provinces in China and directed to Africa to estimate the risk of importation per country. We determined the country's capacity to detect and respond to cases with two indicators: preparedness, using the WHO International Health Regulations Monitoring and Evaluation Framework; and vulnerability, using the Infectious Disease Vulnerability Index. Countries were clustered according to the Chinese regions contributing most to their risk. FINDINGS: Countries with the highest importation risk (ie, Egypt, Algeria, and South Africa) have moderate to high capacity to respond to outbreaks. Countries at moderate risk (ie, Nigeria, Ethiopia, Sudan, Angola, Tanzania, Ghana, and Kenya) have variable capacity and high vulnerability. We identified three clusters of countries that share the same exposure to the risk originating from the provinces of Guangdong, Fujian, and the city of Beijing, respectively. INTERPRETATION: Many countries in Africa are stepping up their preparedness to detect and cope with COVID-19 importations. Resources, intensified surveillance, and capacity building should be urgently prioritised in countries with moderate risk that might be ill-prepared to detect imported cases and to limit onward transmission. FUNDING: EU Framework Programme for Research and Innovation Horizon 2020, Agence Nationale de la Recherche.

An international observational study to assess the impact of the Omicron variant emergence on the clinical epidemiology of COVID-19 in hospitalised patients.

Background: Whilst timely clinical characterisation of infections caused by novel SARS-CoV-2 variants is necessary for evidence-based policy response, individual-level data on infecting variants are typically only available for a minority of patients and settings. Methods: Here, we propose an innovative approach to study changes in COVID-19 hospital presentation and outcomes after the Omicron variant emergence using publicly available population-level data on variant relative frequency to infer SARS-CoV-2 variants likely responsible for clinical cases. We apply this method to data collected by a large international clinical consortium before and after the emergence of the Omicron variant in different countries. Results: Our analysis, that includes more than 100,000 patients from 28 countries, suggests that in many settings patients hospitalised with Omicron variant infection less often presented with commonly reported symptoms compared to patients infected with pre-Omicron variants. Patients with COVID-19 admitted to hospital after Omicron variant emergence had lower mortality compared to patients admitted during the period when Omicron variant was responsible for only a minority of infections (odds ratio in a mixed-effects logistic regression adjusted for likely confounders, 0.67 [95% confidence interval 0.61-0.75]). Qualitatively similar findings were observed in sensitivity analyses with different assumptions on population-level Omicron variant relative frequencies, and in analyses using available individual-level data on infecting variant for a subset of the study population. Conclusions: Although clinical studies with matching viral genomic information should remain a priority, our approach combining publicly available data on variant frequency and a multi-country clinical characterisation dataset with more than 100,000 records allowed analysis of data from a wide range of settings and novel insights on real-world heterogeneity of COVID-19 presentation and clinical outcome. Funding: Bronner P. Gonçalves, Peter Horby, Gail Carson, Piero L. Olliaro, Valeria Balan, Barbara Wanjiru Citarella, and research costs were supported by the UK Foreign, Commonwealth and Development Office (FCDO) and Wellcome [215091/Z/18/Z, 222410/Z/21/Z, 225288/Z/22/Z]; and Janice Caoili and Madiha Hashmi were supported by the UK FCDO and Wellcome [222048/Z/20/Z]. Peter Horby, Gail Carson, Piero L. Olliaro, Kalynn Kennon and Joaquin Baruch were supported by the Bill & Melinda Gates Foundation [OPP1209135]; Laura Merson was supported by University of Oxford's COVID-19 Research Response Fund - with thanks to its donors for their philanthropic support. Matthew Hall was supported by a Li Ka Shing Foundation award to Christophe Fraser. Moritz U.G. Kraemer was supported by the Branco Weiss Fellowship, Google.org, the Oxford Martin School, the Rockefeller Foundation, and the European Union Horizon 2020 project MOOD (#874850). The contents of this publication are the sole responsibility of the authors and do not necessarily reflect the views of the European Commission. Contributions from Srinivas Murthy, Asgar Rishu, Rob Fowler, James Joshua Douglas, François Martin Carrier were supported by CIHR Coronavirus Rapid Research Funding Opportunity OV2170359 and coordinated out of Sunnybrook Research Institute. Contributions from Evert-Jan Wils and David S.Y. Ong were supported by a grant from foundation Bevordering Onderzoek Franciscus; and Andrea Angheben by the Italian Ministry of Health "Fondi Ricerca corrente-L1P6" to IRCCS Ospedale Sacro Cuore-Don Calabria. The data contributions of J.Kenneth Baillie, Malcolm G. Semple, and Ewen M. Harrison were supported by grants from the National Institute for Health Research (NIHR; award CO-CIN-01), the Medical Research Council (MRC; grant MC_PC_19059), and by the NIHR Health Protection Research Unit (HPRU) in Emerging and Zoonotic Infections at University of Liverpool in partnership with Public Health England (PHE) (award 200907), NIHR HPRU in Respiratory Infections at Imperial College London with PHE (award 200927), Liverpool Experimental Cancer Medicine Centre (grant C18616/A25153), NIHR Biomedical Research Centre at Imperial College London (award IS-BRC-1215-20013), and NIHR Clinical Research Network providing infrastructure support. All funders of the ISARIC Clinical Characterisation Group are listed in the appendix.

Impact of spatiotemporal heterogeneity in COVID-19 disease surveillance on epidemiological parameters and case growth rates.

SARS-CoV-2 case data are primary sources for estimating epidemiological parameters and for modelling the dynamics of outbreaks. Understanding biases within case-based data sources used in epidemiological analyses is important as they can detract from the value of these rich datasets. This raises questions of how variations in surveillance can affect the estimation of epidemiological parameters such as the case growth rates. We use standardised line list data of COVID-19 from Argentina, Brazil, Mexico and Colombia to estimate delay distributions of symptom-onset-to-confirmation, -hospitalisation and -death as well as hospitalisation-to-death at high spatial resolutions and throughout time. Using these estimates, we model the biases introduced by the delay from symptom-onset-to-confirmation on national and state level case growth rates (rt) using an adaptation of the Richardson-Lucy deconvolution algorithm. We find significant heterogeneities in the estimation of delay distributions through time and space with delay difference of up to 19 days between epochs at the state level. Further, we find that by changing the spatial scale, estimates of case growth rate can vary by up to 0.13 d-1. Lastly, we find that states with a high variance and/or mean delay in symptom-onset-to-diagnosis also have the largest difference between the rt estimated from raw and deconvolved case counts at the state level. We highlight the importance of high-resolution case-based data in understanding biases in disease reporting and how these biases can be avoided by adjusting case numbers based on empirical delay distributions. Code and openly accessible data to reproduce analyses presented here are available.

Context-specific emergence and growth of the SARS-CoV-2 Delta variant.

The SARS-CoV-2 Delta (Pango lineage B.1.617.2) variant of concern spread globally, causing resurgences of COVID-19 worldwide1,2. The emergence of the Delta variant in the UK occurred on the background of a heterogeneous landscape of immunity and relaxation of non-pharmaceutical interventions. Here we analyse 52,992 SARS-CoV-2 genomes from England together with 93,649 genomes from the rest of the world to reconstruct the emergence of Delta and quantify its introduction to and regional dissemination across England in the context of changing travel and social restrictions. Using analysis of human movement, contact tracing and virus genomic data, we find that the geographic focus of the expansion of Delta shifted from India to a more global pattern in early May 2021. In England, Delta lineages were introduced more than 1,000 times and spread nationally as non-pharmaceutical interventions were relaxed. We find that hotel quarantine for travellers reduced onward transmission from importations; however, the transmission chains that later dominated the Delta wave in England were seeded before travel restrictions were introduced. Increasing inter-regional travel within England drove the nationwide dissemination of Delta, with some cities receiving more than 2,000 observable lineage introductions from elsewhere. Subsequently, increased levels of local population mixing-and not the number of importations-were associated with the faster relative spread of Delta. The invasion dynamics of Delta depended on spatial heterogeneity in contact patterns, and our findings will inform optimal spatial interventions to reduce the transmission of current and future variants of concern, such as Omicron (Pango lineage B.1.1.529).

The Huanan Seafood Wholesale Market in Wuhan was the early epicenter of the COVID-19 pandemic.

Understanding how severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) emerged in 2019 is critical to preventing future zoonotic outbreaks before they become the next pandemic. The Huanan Seafood Wholesale Market in Wuhan, China, was identified as a likely source of cases in early reports, but later this conclusion became controversial. We show here that the earliest known COVID-19 cases from December 2019, including those without reported direct links, were geographically centered on this market. We report that live SARS-CoV-2-susceptible mammals were sold at the market in late 2019 and that within the market, SARS-CoV-2-positive environmental samples were spatially associated with vendors selling live mammals. Although there is insufficient evidence to define upstream events, and exact circumstances remain obscure, our analyses indicate that the emergence of SARS-CoV-2 occurred through the live wildlife trade in China and show that the Huanan market was the epicenter of the COVID-19 pandemic.

Modelling COVID-19.

As the COVID-19 pandemic continues, mathematical epidemiologists share their views on what models reveal about how the disease has spread, the current state of play and what work still needs to be done.

Rapid epidemic expansion of the SARS-CoV-2 Omicron variant in southern Africa.

The SARS-CoV-2 epidemic in southern Africa has been characterized by three distinct waves. The first was associated with a mix of SARS-CoV-2 lineages, while the second and third waves were driven by the Beta (B.1.351) and Delta (B.1.617.2) variants, respectively1-3. In November 2021, genomic surveillance teams in South Africa and Botswana detected a new SARS-CoV-2 variant associated with a rapid resurgence of infections in Gauteng province, South Africa. Within three days of the first genome being uploaded, it was designated a variant of concern (Omicron, B.1.1.529) by the World Health Organization and, within three weeks, had been identified in 87 countries. The Omicron variant is exceptional for carrying over 30 mutations in the spike glycoprotein, which are predicted to influence antibody neutralization and spike function4. Here we describe the genomic profile and early transmission dynamics of Omicron, highlighting the rapid spread in regions with high levels of population immunity.

Recommended reporting items for epidemic forecasting and prediction research: The EPIFORGE 2020 guidelines.

BACKGROUND: The importance of infectious disease epidemic forecasting and prediction research is underscored by decades of communicable disease outbreaks, including COVID-19. Unlike other fields of medical research, such as clinical trials and systematic reviews, no reporting guidelines exist for reporting epidemic forecasting and prediction research despite their utility. We therefore developed the EPIFORGE checklist, a guideline for standardized reporting of epidemic forecasting research. METHODS AND FINDINGS: We developed this checklist using a best-practice process for development of reporting guidelines, involving a Delphi process and broad consultation with an international panel of infectious disease modelers and model end users. The objectives of these guidelines are to improve the consistency, reproducibility, comparability, and quality of epidemic forecasting reporting. The guidelines are not designed to advise scientists on how to perform epidemic forecasting and prediction research, but rather to serve as a standard for reporting critical methodological details of such studies. CONCLUSIONS: These guidelines have been submitted to the EQUATOR network, in addition to hosting by other dedicated webpages to facilitate feedback and journal endorsement.

Genomic epidemiology of SARS-CoV-2 transmission lineages in Ecuador.

Characterisation of severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) genetic diversity through space and time can reveal trends in virus importation and domestic circulation and permit the exploration of questions regarding the early transmission dynamics. Here, we present a detailed description of SARS-CoV-2 genomic epidemiology in Ecuador, one of the hardest hit countries during the early stages of the coronavirus-19 pandemic. We generated and analysed 160 whole genome sequences sampled from all provinces of Ecuador in 2020. Molecular clock and phylogeographic analysis of these sequences in the context of global SARS-CoV-2 diversity enable us to identify and characterise individual transmission lineages within Ecuador, explore their spatiotemporal distributions, and consider their introduction and domestic circulation. Our results reveal a pattern of multiple international importations across the country, with apparent differences between key provinces. Transmission lineages were mostly introduced before the implementation of non-pharmaceutical interventions, with differential degrees of persistence and national dissemination.

Assessing the impact of COVID-19 border restrictions on dengue transmission in Yunnan Province, China: an observational epidemiological and phylogenetic analysis.

BACKGROUND: In response to the COVID-19 pandemic, China implemented strict restrictions on cross-border travel to prevent disease importation. Yunnan, a Chinese province that borders dengue-endemic countries in Southeast Asia, experienced unprecedented reduction in dengue, from 6840 recorded cases in 2019 to 260 in 2020. METHODS: Using a combination of epidemiological and virus genomic data, collected from 2013 to 2020 in Yunnan and neighbouring countries, we conduct a series of analyses to characterise the role of virus importation in driving dengue dynamics in Yunnan and assess the association between recent international travel restrictions and the decline in dengue reported in Yunnan in 2020. FINDINGS: We find strong evidence that dengue incidence between 2013-2019 in Yunnan was closely linked with international importation of cases. A 0-2 month lag in incidence not explained by seasonal differences, absence of local transmission in the winter, effective reproductive numbers < 1 (as estimated independently using genetic data) and diverse cosmopolitan dengue virus phylogenies all suggest dengue is non-endemic in Yunnan. Using a multivariate statistical model we show that the substantial decline in dengue incidence observed in Yunnan in 2020 but not in neighbouring countries is closely associated with the timing of international travel restrictions, even after accounting for other environmental drivers of dengue incidence. INTERPRETATION: We conclude that Yunnan is a regional sink for DENV lineage movement and that border restrictions may have substantially reduced dengue burden in 2020, potentially averting thousands of cases. Targeted testing and surveillance of travelers returning from high-risk areas could help to inform public health strategies to minimise or even eliminate dengue outbreaks in non-endemic settings like southern China. FUNDING: Funding for this study was provided by National Key Research and Development Program of China, Beijing Science and Technology Planning Project (Z201100005420010); Beijing Natural Science Foundation (JQ18025); Beijing Advanced Innovation Program for Land Surface Science; National Natural Science Foundation of China (82073616); Young Elite Scientist Sponsorship Program by CAST (YESS) (2018QNRC001); H.T., O.P.G. and M.U.G.K. acknowledge support from the Oxford Martin School. O.J.B was supported by a Wellcome Trust Sir Henry Wellcome Fellowship (206471/Z/17/Z). Chinese translation of the abstract (Appendix 2).

Pneumonia of unknown aetiology in Wuhan, China: potential for international spread via commercial air travel.

There is currently an outbreak of pneumonia of unknown aetiology in Wuhan, China. Although there are still several unanswered questions about this infection, we evaluate the potential for international dissemination of this disease via commercial air travel should the outbreak continue.

Spatiotemporal invasion dynamics of SARS-CoV-2 lineage B.1.1.7 emergence.

Understanding the causes and consequences of the emergence of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) variants of concern is crucial to pandemic control yet difficult to achieve because they arise in the context of variable human behavior and immunity. We investigated the spatial invasion dynamics of lineage B.1.1.7 by jointly analyzing UK human mobility, virus genomes, and community-based polymerase chain reaction data. We identified a multistage spatial invasion process in which early B.1.1.7 growth rates were associated with mobility and asymmetric lineage export from a dominant source location, enhancing the effects of B.1.1.7's increased intrinsic transmissibility. We further explored how B.1.1.7 spread was shaped by nonpharmaceutical interventions and spatial variation in previous attack rates. Our findings show that careful accounting of the behavioral and epidemiological context within which variants of concern emerge is necessary to interpret correctly their observed relative growth rates.

Transmission of SARS-CoV-2 before and after symptom onset: impact of nonpharmaceutical interventions in China.

Nonpharmaceutical interventions, such as contact tracing and quarantine, have been the primary means of controlling the spread of SARS-CoV-2; however, it remains uncertain which interventions are most effective at reducing transmission at the population level. Using serial interval data from before and after the rollout of nonpharmaceutical interventions in China, we estimate that the relative frequency of presymptomatic transmission increased from 34% before the rollout to 71% afterward. The shift toward earlier transmission indicates a disproportionate reduction in transmission post-symptom onset. We estimate that, following the rollout of nonpharmaceutical interventions, transmission post-symptom onset was reduced by 82% whereas presymptomatic transmission decreased by only 16%. The observation that only one-third of transmission was presymptomatic at baseline, combined with the finding that NPIs reduced presymptomatic transmission by less than 20%, suggests that the overall impact of NPIs was driven in large part by reductions in transmission following symptom onset. This implies that interventions which limit opportunities for transmission in the later stages of infection, such as contact tracing and isolation, are particularly important for control of SARS-CoV-2. Interventions which specifically reduce opportunities for presymptomatic transmission, such as quarantine of asymptomatic contacts, are likely to have smaller, but non-negligible, effects on overall transmission.

Genomics and epidemiology of the P.1 SARS-CoV-2 lineage in Manaus, Brazil.

Cases of severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) infection in Manaus, Brazil, resurged in late 2020 despite previously high levels of infection. Genome sequencing of viruses sampled in Manaus between November 2020 and January 2021 revealed the emergence and circulation of a novel SARS-CoV-2 variant of concern. Lineage P.1 acquired 17 mutations, including a trio in the spike protein (K417T, E484K, and N501Y) associated with increased binding to the human ACE2 (angiotensin-converting enzyme 2) receptor. Molecular clock analysis shows that P.1 emergence occurred around mid-November 2020 and was preceded by a period of faster molecular evolution. Using a two-category dynamical model that integrates genomic and mortality data, we estimate that P.1 may be 1.7- to 2.4-fold more transmissible and that previous (non-P.1) infection provides 54 to 79% of the protection against infection with P.1 that it provides against non-P.1 lineages. Enhanced global genomic surveillance of variants of concern, which may exhibit increased transmissibility and/or immune evasion, is critical to accelerate pandemic responsiveness.